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ABSTRACT 

Aims. The main goal of this study is to detect the stellar overdensity associated with the Perseus arm in the anticenter direction. 
Methods. We used the physical parameters derived from Stromgren photometric data to compute the surface density distribution 
as a function of galactocentric distance for different samples of intermediate young stars. The radial distribution of the interstellar 
absorption has also been derived. 

Results. We detected the Perseus arm stellar overdensity at 1.6±0.2 kpc from the Sun with a significance of 4-5cr and a surface density 
amplitude of around 10%, slightly depending on the sample used. Values for the radial scale length of the Galactic disk have been 
simultaneously fitted obtaining values in the range [2.9,3.5] kpc for the population of the B4-A1 stars. Moreover, the interstellar visual 
absorption distribution is congruent with a dust layer in front of the Perseus arm. 

Conclusions. This is the first time that the presence of the Perseus arm stellar overdensity has been detected through individual star 
counts, and its location matches a variation in the dust distribution. The offset between the dust lane and the overdensity indicates that 
the Perseus arm is placed inside the co-rotation radius of the Milky Way spiral pattern. 
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1. Introduction 


The spiral arm structure in the Galactic disk is an important 
component when studying the morphology and dynamics of the 
Milky Way. However, we still lack a complete picture that de¬ 
scribes the nature, origin, and evol ution of this s tructure. Are 
there two (iDrimmefeOOOt) or four dRusseill 120031) spiral arms? 
Are their main component s g as or stars? Maybe, as was re¬ 
cently suggested bv IBeni aminl 1 2008 ), there are two major spiral 
arms (Scutum-Centaurus and Perseus) with higher stellar densi¬ 
ties and two minor arms (Sagittarius and Norma) mainly filled 
with gas and star forming regions. It is also important to under¬ 
stand how the stars interact with the arms: do they move through 
the arms (density wave me cha nism, Lin & S hulll964l) . along the 
arms (manifold s. iRomero-Gomez et al.l 20071) . or with the arms 
(material arms, [Grand et al.ll2012l) ? 

The Perseus spiral arm, with its outer structure placed near 
the Sun, is an excellent platform from which to undertake a 
study of its nature. It has been studied considering different 
tracers such as HI neutral gas (lLindbladlll967l) : large-scale CO 
distribution surveys (iDame et all 2001 ): a compilation of open 
clusters and associati ons and CO surveys of m olecular clouds 
(VazcjuezeUil. 2008 ); star forming complexes (iRussei IHqqI 
iFoster & Brunt 2014 ); or masers associate d with young, high- 
mass stars ( Xu et al . 2006HReid et al.l2014h . Some of these stud¬ 


ies trace the arm only in the second quadrant, and others trace it 
in the third quadrant, but few analyses link both. Furthermore, 
few of them have enough data pointing towards £ = 180°, that is 
towards the anticenter. Recently, [Reid et akl (120141) published a 
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substantial compilation of over 100 high-precision trigonomet¬ 
ric distances to several masers using VLBI observations, in other 
words tracing the location of star forming regions. Around 20 
masers are located in the Perseus spiral arm, with five of them 
very close to the anticenter. Their model for the spiral arms lo¬ 
cates the Perseus arm at 2 kpc at around £ = 180°, with the five 
masers placed at slightly smaller galactocentric radius than the 
fit. IValleel (120141 ) published a master catalog of the observed 
tangents to the Galaxy’s spiral arms in order to fit a four-arm 
model. Using different arm tracers, he located the “Perseus ori¬ 
gin arm” near the Galactic center at £ — 338°. The fitted model 
allows us to extrapolate the position of the Perseus arm at about 
2 kpc in the anticenter direction. Very interesting is his attempt 
to quantify the offset between different tracers, which favors the 
interpretation of the data in terms of the density wave theory. 

The present study aims to trace the Perseus arm by study¬ 
ing both the radial stellar density variation toward the anticenter 
with simple star counts, and by deriving the distribution of the ra¬ 
dial velocity components through this direction. Both arm struc¬ 
ture and kinematics are essential issues in our study. To trace 
the arm we use intermediate young stars with effective tempera¬ 
tures in the range [15000,9000] K (B4-A1 spectral type). These 
stars are excellent tracers of this overdensity as 1) they are bright 
enough to reach large distances from the Sun, and 2) they are 
old enough to have had time to respond to the spiral arm po¬ 
tential perturbation. We explicitly avoid very young 0-B3 and 
Cepheids since we aim to determine the location of the mass 
peak (i.e., the potential minimum) associated with the Perseus 
arm and not the peak of star formation. In the density wave sce¬ 
nario, one would expect a sequence (in distance or angle) starting 
with dust and then star formation, but offset with respect to the 
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potential dRobe 

rts|[i969 

). This scenario is in agreement with the 

recent work of 

Vallee 

2014j) where the hot dust (with masers 


and newborn stars) seems to peak near the inner arm edge while 
the stars are all over the arms. Furthermore, to analyze the inter¬ 
action between the arm and the stellar component, the selected 
B4-A1 stellar tracers are young enough so their intrinsic veloc¬ 
ity is still small, so that their response to a perturbation would 
be stronger and therefore easier to detect. Our tracers -late B- or 
early A-type stars- are expected to show a density variation due 
to the presence of a perturbation. As mentioned, the most suit¬ 
able direction to undertake this study on Galactic structure and 
kinematics is towards the anticenter and the Perseus arm. First, 
this direction presents lower interstellar extinction than the di¬ 
rection pointing to the Galactic center. Second, although slightly 
depending on the pitch angle, this direction is the one expected 
to have the smallest distance to the Perseus arm stellar compo¬ 
nent. In addition, and more importantly, this direction was se¬ 
lected because the Galactic rotation would be negligible in the 
stellar radial velocity component spectroscopically derived for 
this young population, and it would directly reflect the velocity 
perturbation introduced by the spiral arm. 

In the present paper we present the analysis of the stellar sur¬ 
face density through accurate distances and interstellar absorp¬ 
tion derived from Stromgren photometry. In a second paper of 
this series, work on the spiral arm kinematic perturbation will 
be provided utilizing stellar radial velocity for a subsample of 
selected young stars. 

In Sect. [2] we overview the main characterist ics o f our pho- 
tometric surv ey published and characterized in IMonguio et alJ 
(I2013L |20I4|) . In Sect. [2 we carefully select the distance lim¬ 
ited samples required for this study. All observational biases and 
constraints are evaluated and accounted for. This analysis allows 
us to estimate in Sect. [4] the local peaks in the surface density 
distribution compared to a pure exponential profile. In Sect. 0 
we characterize the stellar overdensity associated with Perseus, 
its distance from the Sun, the spiral arm amplitude, and its sig¬ 
nificance. The differential interstellar absorption distribution and 
its relation with the Perseus spiral arm are discussed in Sect. [6] 
Finally, in Sect. [7J we summarize the main results and conclu¬ 
sions. 

2. Data 

As published in iMonguio et alJ 1 2013) . the authors carried out 
a uvbyH/3 Stromgren photometric survey covering 16° 2 in the 
anticenter direction using the Wide Field Camera at the Isaac 
Newton Telescope (INT, La Palma). This is the natural photo¬ 
metric system for identifying young stars and obtaining accurate 
estimates of individual distances and ages. The survey was cen¬ 
tered slightly below the plane in order to take into account the 
warp, and covers Galactic longitudes from C ~177° to t ~ 183° 
and Galactic latitudes from b -2° to b ~ L5. The calibra¬ 

tion to the standard system was undertaken using open clusters. 
We created a main catalog of 35974 stars with all Stromgren 
indexes and a more extended one with 96980 stars with partial 
data. The inner 8° 2 reach ~90% completeness at V ~ 17 m , while 
the outer sky area of ~8° 2 , mostly observed with only one point¬ 
ing, reaches this completeness at V ~ 15™5. Photometric inter¬ 
nal precisions around 0.01-0.02™ for stars brighter than V — 16 m 
were obtained, increasing to O'."05 for some indexes and fainter 
stars ( V =1 8 -19 m ). _ 

In iMonguio et al. ( 201 4) we describe in detail two different 
approaches implemented to compute the stellar physical parame¬ 
ters (SPP) for these stars. The first uses available pre-Hipparcos 


photometric calibrations from the 1980s (empirical calibration 
method, EC) based on previous cluster data and trigonometric 
parallaxes. This procedure follows two steps: 1) the classifica¬ 
tion of the stars in different photometric regions using extinction- 
free indexes ([ci], [mi], Hj3, [u - b]), and 2) for each region, the 
interpolation in empirical calibration sequences to obtain the in¬ 
trinsic indexes and the absolute magnitudes from which inter¬ 
stellar extinction and distances can be computed. The second 
method, which we developed, is based on the most accurate at¬ 
mospheric models and evolutionary tracks available at present 
(model based method, MB). It starts with a 3D fit using three 
extinction free indexes ([ci], [mi], ///)), and their photometric er¬ 
rors. Then, the interpolation in evolutionary tracks provides the 
remaining physical parameters. Both methods were optimized 
for hot stars using the comparison with Hipparcos data. In both 
cases, individual errors on SPP were computed through Monte 
Carlo random realizations. 

The MB method provided distance. My, Ay, ( b - v)o, T e ff, 
and log#, among other SPP for stars with T e ff >7000K. On 
the other hand, the EC method provided distance. My, Ay, and 
(b - y)o for stars with spectral types in the range B0-A9. By 
comparing these two sets of SPP data, a clear discrepancy of 
about ~20% in distance was detected and discussed. The catalog, 
published in IMonguio et al] (120141) . also provided quality flags, 
some of them related to the phot ometric c lassification required 
by the classical EC method (IStrdmgrenl 19661) to disentangle be¬ 
tween early- and late-type stars. It was evident that the separation 
in the [ci] - [mi] space is well defined and works properly for 
stars with low photometric errors, but the gap between regions 
blurs when faint stars with larger photometric errors are consid¬ 
ered. The IPHAS and 2MASS data, when available, were used 
to flag emission line stars and to check the coherence of some 
physical parameters derived using visual and infrared data, re¬ 
spectively. The fraction of binary stars in ou r samp l e may b e im- 
portant, in particular for hot stars. Thus in lMonguio et al.l (120141) 
the effects of the secondary on the final accuracy of the photo¬ 
metric indexes and distance were also discussed. 

3. Distance complete samples 

The available data -full uvbyH/3 photometry for 35974 stars in 
the anticenter direction- were cleaned and organized in different 
working samples by carefully evaluating both the quality of the 
data and the apparent magnitude and distance completeness (see 
Fig-EQ- In a fi rst st ep, the ava ilab le IPHAS initial data release 
(IDR) dGonzalez-Solares et al .1120081) . combined with our H[3 in¬ 
dex, allowed us to detect and reject 1073 emission line stars. 
Then each approach to compute the physical parameters (i.e., 
MB and EC) required its own strategy. We created two differ¬ 
ent working samples: the first, hereafter MB-S1, with the 12957 
stars with T e ff > 7000K and SPP from the MB method, and 
the second, hereafter EC-S1, with 11524 B0-A9 stars, with SPP 
derived using the EC method. Following the MB strategy, we 
considered as outliers the stars that lie 5<x outside the MB at¬ 
mospheric grid during the 3D fit, taking into account the photo¬ 
metric errors in [ci] - [mi] - Hji. We also removed stars with 
doubtful classification between photometric regions using N S id e 
and N reg flags. These flags, N reg for the EC method and /V v „/ ( . for 
the MB method, indicate the probability for a star to belong to 
earlier or later photometric regionsand were computed through 
Monte Carlo simulations (see iMonguio et akll2014l for details). 
Coherence with 2MASS data allowed us to reject some stars that 
cannot be allocated reliably to either early (B0-B9) or late (A3- 
A9) types. The two resulting “clean” samples contain 9781 (for 
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Fig. 1: Procedure for the generation of the working samples. 
From top to bottom, the number of stars in the samples decreases 
while the quality of the individual physical parameters increases. 
The left side shows the samples according to the MB method, 
and the right side with EC. 


MB-S2) and 9622 stars (for EC-S2). A new subset was built con¬ 
taining the stars that simultaneously belong to the young popula¬ 
tion following both methods. This was done by checking the co¬ 
herence between the two sets of physical parameters. This clean¬ 
est sample contains 8328 stars and is named CS-MB or CS-EC 
depending on whether their physical parameters were computed 
using the MB or EC method (see Fig.Q}. This sample contains 
stars with more accurate SPP data, but has fewer statistics owing 
to the lack of stars rejected during the cleaning process. From 
now on we will distinguish between the subsamples containing 
the stars in the inner sky area (in suffix), from those where the 
stars located in the outer sky area are added (all suffix). Whereas 
the first subsample reaches fainter limiting magnitude, the sec¬ 
ond contains more stars but is complete up to a brighter limiting 
magnitude. In Table |T] we give the number of stars for each of 
the subsamples. 

Since we wanted to derive how the stellar density varies with 
the galactocentric distance, we needed to ensure that the features 
observed are due to a physical and real overdensity and not to 
observational or selection biases. To account for this, we created 
distance complete samples by selecting ranges of absolute mag¬ 
nitude, that ensured us distance completeness up to a given helio¬ 
centric distance. For each sample we estimated 1) the maximum 
visual limiting magnitude Vn m , 2) the magnitude at which some 
stars were saturated V sat , and 3) the maximum and minimum vi¬ 
sual interstellar absorption at a given distance. The Vu m for each 


of the samples, assuming a 90% completeness, was estimated by 
computing the mean of the magnitudes at the peak star counts in 
a magnitude histogram and its two adjacent bins, before and af¬ 
ter the peak, weighted by the number of stars in each bin. In 
[Monguio et al.l (1201 3l) we checked that this estimate is adequate 
by comparing the V magnitude distribution of the stars with all 
the available photometric indexes, and with the distribution ob¬ 
served for stars measured in V that reached a significant fainter 
limiting magnitude. To ensure an adequate range before and after 
the Perseus arm in the anticenter direction, we selected samples 
complete up to 3 kpc from the Sun. The maximum interstellar 
absorption at 3 kpc was estimated from A Vmax (3 kpc)=Ay + cta v ■ 
This value was then used to compute the limit in the intrinsic star 
brightness: M VUm ( 3 kpc) = V Um - A Vmax (3kpc) - 5 log(3 kpc) + 5. 
In addition, stars brighter than V„, = ll m may suffer from 
saturation problems ([Monguio et al.ll2013l) . Thus, we estimated 
and set a completeness limit in the absolute brightness -named 
Myu m (\ ,2kpc)~ for the closest distance, imposed to be 1.2kpc, 
using both V sat and the interstellar absorption at this distance. 
Following this strategy, that is, by selecting the stars with the 
absolute magnitude inside the computed range, all the resulting 
samples were converted to distance complete samples in the dis¬ 
tance range between 1.2 and 3 kpc. Table[j]shows all the values 
used to generate each sample together with the number of stars 
included in each of them. Figure [2] shows an example of the ap¬ 
parent visual magnitude range needed to cover the distance range 
[1.2,3] kpc for a star with a fixed absolute magnitude. As an ex¬ 
ample, this figure shows the limits on apparent visual magnitude 
assumed for the MB-S1, MB-S2, and CS-MB samples. We can 
see how these samples have complete populations inside the he¬ 
liocentric distance range selected. 

Finally, we checked that the distribution of effective temper¬ 
ature and spectral type included in each of our samples do not 
change as a function of distance within the [1.2,3] kpc range. 
Several two-sample Kolmogorov-Smirnov (KS) tests were done 
comparing the luminosity function at different distances with the 
distribution at the solar neighborhood (obtained in Sect. 14.21 see 
Fig-0- In a H the cases, the KS tests are consistent with the as¬ 
sumption that the underlying luminosity functions arise from the 
same distribution (see Table 0. We observed that only for the 
farthest subsamples of MB-S1 and MB-S2 we do have small val¬ 
ues of the p-values, but not enough to reject the null hypothesis 
at 0.01 level. 

Booking at Table Q] we see that the samples considering all 
the anticenter observed sky area, as the outer area has a brighter 
limiting magnitude, have a very narrow range in absolute magni¬ 
tude, which means that the final number of stars in those samples 
is not high enough to undertake the study of the stellar overden¬ 
sity induced by the Perseus arm with good statistical significance 
(Sect.0. 

4. Radial density distribution 

The surface density was chosen over volume density to be the 
best parameter with which to trace the stellar distribution con¬ 
sidering the information available, that is, a set of distance com¬ 
plete samples inside a solid angle of 8° 2 or 16° 2 in the anticen¬ 
ter direction (Sect. 0. This is a better choice since it gives an 
estimate of the total disk mass at a given distance, also taking 
into account the change in the mid-plane due to the warp, as¬ 
suming that the disk scale height varies slowly with distance. 
All the factors taken into account in the derivation of this pa¬ 
rameter from the data are described in Sect. 14.11 The radial bin 
size used to derive the galactocentric variation of this parameter 
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Table 1: Ay, V/ lm , and My/,,,, for all samples at the bright (min) and faint (max) ends. Spectral type according to the Mvn m are shown. 
N is the number of stars for the absolute magnitude limited samples, i.e., Myn m ( 1.2kpc) < My < Myu m ( 3kpc), also indicating 
which of them are between 1.2 and 3 kpc. N/N^i is the ratio between the number stars (between 1.2 and 3 kpc) of each sample and 
the SI sample, always using the My range of the sample. 



MB 

EC 


in 

all 

in 

all 


SI I S2 I CS 

SI | S2 | CS 

SI I S2 I CS 

SI | S2 | CS 

\/min 

V lim 

11 

\/max 

V lim 

16.7 | 16.7 | 16.3 | 15.7 | 15.7 | 15.3 

16.7 | 16.7 | 16.3 | 15.7 | 15.7 | 15.3 

A Vmin ( 1 • 2kpC ) 

1.5 

1.8 

Av, !MJ (3kpc) 

3.1 

3.2 

M v/i „(1.2kpc) 

-0.9 

-1.2 

SP,„,„ 

B4 

B3.5 

Mvumi 3kpc) 

1.2 

1.2 

0.8 

0.2 

0.2 

-0.2 

1.1 

1.1 

0.7 

0.1 

0.1 

-0.3 

SP 

VJJ - max 

Al 

Al 

A0 

B8.5 

B8.5 

B7.5 

A0 

A0 

B9.5 

B8.5 

B8.5 

B7 

N to , 

7763 

6061 

4998 

12957 

9781 

8101 

6756 

5943 

4998 

11524 

9622 

8101 

N (My/„„) 

2209 

1565 

672 

1587 

935 

278 

821 

653 

361 

446 

289 

133 

N(1.2-3kpc) 

851 

733 

349 

403 

326 

106 

441 

378 

185 

151 

111 

33 

N/Nsi 

1.00 

0.86 

0.65 

1.00 

0.81 

0.63 

1.00 

0.86 

0.82 

1.00 

0.74 

0.57 


Table 2: Two sample KS test results of the comparison between solar luminosity function and luminosity function at different 
distances ranges. We provide the KS statistic D nun2 , the number of stars for each of the samples (n \ for the solar sample, hi for the 
sample at each distance), and the p-value. 




1.2-1.5 kpc 

Dn i ,«2 ^2 P 

1.5-2.0 kpc 

Dfi 1 ,«2 ^2 

P 

2.0-2.5 kpc 

Dn [ ,«2 ^2 

P 

2.5-3.0kpc 

D n j ,n 2 ^2 

P 

MB-Slin 

343 

0.10 

91 

0.44 

0.09 

231 

0.20 

0.06 

236 

0.73 

0.12 

293 

0.03 

MB-S2in 

343 

0.11 

81 

0.44 

0.09 

200 

0.27 

0.08 

202 

0.45 

0.12 

250 

0.03 

CS-MBin 

195 

0.16 

40 

0.37 

0.09 

89 

0.75 

0.15 

98 

0.09 

0.12 

122 

0.23 

EC-Slin 

36 

0.15 

62 

0.61 

0.16 

129 

0.43 

0.10 

134 

0.90 

0.16 

116 

0.46 

EC-S2in 

36 

0.21 

55 

0.27 

0.16 

116 

0.47 

0.13 

117 

0.69 

0.14 

90 

0.67 

CS-ECin 

18 

0.24 

20 

0.56 

0.25 

62 

0.31 

0.14 

56 

0.94 

0.13 

47 

0.98 



Fig. 2: Observed apparent magnitude for different absolute mag¬ 
nitude stars lo cated at different distances. We have assumed an 
absorption of Ay{r) + <r Av ( r) obtained from our data using the 
MB method. Vertical lines show the established limits at 1.2 
and 3 kpc. Horizontal lines show the Vu m computed for different 
working samples. In red, V / '"" A and AF'”" A forMB-Sl and MB-S2 
samples. In blue, V™" x and M™* forCS-MB. In green, V™" and 
M'ylim com puted using Ay(r). 


was computed using the method proposed bv lKnu thl d2006i) . As 


stated by the author, this method is optimized to find substruc¬ 
ture in the datiQ. The radial scale length of the Galactic disk in 
the anticenter direction and the possible overdensity associated 
with Perseus were obtained by fitting an exponential function 
to the working samples described in Table |T| Initially, the two 
parameters of this function were fitted: the radial scale length 
( Iir ) and the surface density in the solar neighborhood (T 0 ,) of 
the population represented by these working samples. The fit 
highly depends on the zero point of the distribution Z G . Thus, 
to substantially increase both the range in distance covered and 
the significance of the fit, we decided to derive this parameter 
independently using a local sample (see Sect. 14. 2I >. Then these 
values were used in Sect. !4.3l to do the fits, where the departure 
of the exponential surface density distribution and the existence 
of possible maxima are discussed in terms of the location and 
existence of the overdensity associated with the Perseus spiral 
arm. 


4. 1. The computation of the stellar surface density 

To derive the stellar surface density we included the sky area 
surveyed, the effect of the scale height of the underlying thin 
disk population being considered, as well as the presence of the 
Galactic warp. The surface density for each distance bin was 
computed as 


1 Nonetheless, we checked that a simpler procedure such as the con¬ 
stant radial bin of 200 pc also provides congruent results with equivalent 
conclusions. 
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Z(r k ) = — J —, 

Sk tr F ZJ 


(1) 


where is the number of stars in a given radial bin k; Sk is 
the disk projected surface for that radial bin computed from the 
Galactic latitude £ covered, the mean distance to the Sun r/ ( of 
the bin, and the distance width size of the bin in parsecs; Fzj 
is a correction factor for each star that takes into account the 
fact that our observed solid angle does not cover the full vertical 
cylinder. This factor is assumed to depend on the vertical den¬ 
sity distribution, modeled as sech 2 (z/h z ) (Ivan der Kruit & Searle l 
1 1981 )1. with h z being the scale height of the disk. As is known, 
h z depends on the age of the population. Here, as an approxi¬ 
mation, we model it as depending on the visual absolute magni¬ 
tude of the star. Values for h, are not well determined, whereas 
iReedl (120001) gave h 7 =2 5 -65 pc for OB-type stars (mainly O- 
B2), and iMafz-Apellanizl (12001 1) obtai ned /z 7 -34.2±2.5pc for a 
sample of 0-B5 stars. More recently, iKong & Zhul (120081) . us¬ 
ing Hipparcos data, estimated values_around h -- 10 3.1 +3.0pc 
for the full range of OB stars. ICzekai et al.l ( 20141) adopted 
/i.= 130 pc for stars with ages younger than 0.15Gyr whereas 
values up to 260 pc were assumed for stars with ages between 
0.15 and 1 Gyr. These simulations provide age ranges around 
50+20 Myr for B5-type stars (J e ff ~ 15000 K), while age ranges 
of 500+200Myr are obtained for A5-type stars (Jeff ~ 8000 K). 
Taking into account these estimations, we adopted the relation 
/i.(pc) = 36.8 ■ My + 130.9, which is equivalent to assuming 
100pc for B5-type stars and 200 pc for an A5-type star. It must 
be taken into account that a different relation between the scale 
height and the intrinsic brightness would change the zero point 
of the stellar density distribution, and thus the corresponding 
surface density at the Sun’s position S 0 - The Fzj factor also 
takes into account the vertical range ( z m i n ,i, z m ax,i) covered at each 
galactocentric distance. Keeping this in mind, this factor is com¬ 
puted following 


Fzj(r ) = 


£ 


sech 


(- 

U 1 


~ Zw(r) )dz 


J_7 sech2 (l$S) dz ' 


( 2 ) 


where we also considered the position of the warp (zw(r)) . From 


.. rj). 1 

observational data (2MASS and HI data) iMomanv et all ( 2006 ) 
obtained bw —0.5° at £ = 180°, so we can assume that at differ¬ 
ent distances the zw, where the star density of the disk is maxi¬ 
mum, can be computed as zw - r ■ tan b w . 


4.2. Surface density in the solar neighborhood 

We used the Hauck & Mermilliod ( 1998 ?) catalog of Stromgren 
photometry to compute the surface density at the Sun’s posi¬ 
tion for young stars. This process was done following the same 
methodology as we used for the anticenter stars, that is, the 
same computation of SPP, and the same process used to clean 
the sample s (Sects. |2] and [3] >. First we checked the complete¬ 
ness of the lHauck & Mermilliod ( 11998b catalog in terms of the 
visual apparent magnitude. For that, we cross-matched this cata¬ 
log with the Hipparcos catalog (ESA 1997). We verified that the 
completeness is above 95% up to 6™5 for OB stars. Stromgren 
photometric data allowed us to compute their physical param¬ 
eters, selecting only those with T e ff > 7000K for MB, and 
B0-A9 for EC. To mimic the samples in the anticenter we se¬ 
lected the same My ranges, and computed the distance limit 
for which we can ensure completeness (assuming in this case 
Ay = 0): ru m - l0 (y '™"” Mv '"“ + W 5 . We repeated the computations 



Fig. 3: Surface density at the Sun’s position for 0.2 m My bins, 
using both Vn m - 6 m (in red) and V/, m = 6‘."4 (in blue). In black 
we have the same values as obtained from the Besanfon Galaxy 
model. 


using Vum = 6 m and V / /,„, =6 r ."4 and obtained very similar re¬ 
sults. We considered the effect of the warp to be negligible in 
the solar neighborhood. We also checked that different values 
for the distance of the Sun above the Gala ctic plane gave results 
within the err or bars; Zn =1 5 pc fr om iRobin et all d2003h and 
Zo=26pc from M aiaess et al J (l2009t) . We used the same depen¬ 
dence of the scale height function on absolute magnitude as for 
the density in the anticenter direction (Sect. 14.11 ). In this case, 
the limits in z follow the surface of an sphere around the Sun’s 

position, Zmax = Z m in = xj rj im - x 2 - y 2 , and the projected sur¬ 


face density is then easily modeled as Sk = tr r 2 jm . The obtained 
distribution for different magnitudes ranges of (). n 2 is plotted in 
Fig .[3] For comparison, the values resulting fro m m odel B of 
the new Besangon Galaxy Model (ICzekai et al.ll2014l) are over¬ 
plotted. These values were derived fitting the model to the full 
sky Tycho catalogue. In Fig. [3] we observe that the distribution 
is very similar in the whole range of absolute magnitudes, thus 
confirming the robustness of the derivation of the local surface 
density for the young population, i.e., the zero point of our ex¬ 
ponential fit to the anticenter data. When we used the My limits 
established for our samples (see Table 0}. the values for the lo¬ 
cal surface density obtained were (3.25±0.13)-10~ 2 for MB-Slin 
andMB-S2in, (1.70±0.07)-l(T 2 forCS-MBin, (3.49+0.12)-1(T 2 
for EC-Slin and EC-S2in, and (1.41+0.06)-1 O' 2 */pc 2 for CS- 
ECin. Their errors were computed as 


02 


yfii 

Z’UFzjSk/n 


(3) 


4.3. The radial scale length and the detection of the Perseus 
arm 

The radial surface density distribution was computed for all the 
samples and an exponential function was fitted to each of them 
by adopting the zero point values for the local surface den¬ 
sity obtained in the previous section. The radial scale length 
of the Galactic disk in the anticenter direction is presented in 
Fig.[4] The scale lengths obtained for the fits using MB distances 
(Fig. 13 top panel) are in the range !ir —[3.5-4.3] kpc. More 
importantly, we can clearly recognize the overdensity associ¬ 
ated with the Perseus arm at around 1.6 kpc (depending slightly 
on the sample), within the distance range at which the sam¬ 
ples are complete, that is, between 1.2 and 3 kpc from the Sun. 
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The decline after 3 kpc is produced by the limiting magnitude 
of our samples, since we carefully checked that the complete¬ 
ness is only ensured up to this distance. The results for the EC 
method are presented in the bottom panels of Fig. [4] As dis¬ 
cussed in Sect. [2] the photometric distances derived from this 
method are clearly biased, in the sense that EC distances can 
be 20% shorter than the MB distances, the latter being more 
congruent with Hip parcos parallaxes (the detailed comparison 
is given in lMonguio et al.ll2014l) . As a consequence of this bias, 
the overdensity bump associated with the Perseus arm is also 
present, but closer to us and overlapping with the short distance 
limit at 1.2 kpc imposed to avoid saturation effects from bright 
stars. Some attempts were made to fit the points beyond the arm, 
but the results obtained for the radial scale length are unreal¬ 
istic (/2/j=[1.7,1.9]kpc). For the samples using all the sky area, 
the number of stars in the range [1.2-3] kpc was too low, since 
the limiting magnitude was too bright, and the fainter stars in¬ 
cluded in the samples have My-0.2,-0.2, i.e., B7-B8.5 stars (see 
Table Q}. Thus, the features related to the overdensity, although 
present, were less statistically significant than for the in samples. 

The next step was to re-compute the exponential fitting, 
avoiding the points close to the peak detected, that is between 
1.4 and 2.0 kpc (see Fig.0. In this fit, a slight decrease of the ra¬ 
dial scale length until Iir =2.9+0.1/0.2 kpc for both MB-S1 and 
MB-S2 samples is obtained. For the cleanest sample CS-MB we 
obtain a slightly larger value of Iir =3.5+0.5 kpc (right panel of 
Fig. [3- In this case, the uncertainty is significantly larger since 
the working sample contains fewer stars. In addition, the stars 
in this sample have a brighter apparent limiting magnitude with 
respect to the previous ones, so their corresponding My range is 
shifted to intrinsically brighter objects, i.e., it contains a slightly 
younger population on average. This result is congruent with re¬ 
cent data o btaine d from simulations such as those presented by 
iBird et akl d2013l) . showing that the radial scale length increases 
when the population considered is younger. 

The three samples using MB distances show a clear peak 
around 1.6 kpc (Figs. [4] and 0). This overdensity is clearer for 
those samples having larger number of stars (MB-S1 and MB- 
S2, left and central panels). Although they have a more robust 
statistical detection, they also contain values that are less accu¬ 
rate for the SPP and some contamination from adjacent photo¬ 
metric regions could also be present (see Sect. 0. Thus we con¬ 
sider that it is the whole set of fits presented in the top rows of 
Figs. [4] and [3 using the three samples, that allows us to con¬ 
firm the detection of the stellar overdensity associated with the 
Perseus spiral arm. 


5. The Perseus arm overdensity 

We performed^ 2 tests for the different fits presented in the pre¬ 
vious sections to quantitatively evaluate the significance of the 
detected Perseus arm overdensity. The number of stars observed 
{n° k bs , with k = l, m) for each k distance bin between 1.2 and 
3 kpc was compared with the distribution coming from the re¬ 
sulting exponential fit (see Fig. [3- To compare the model and 
the observations we use as a first approximation the number of 
corresponding stars at the central distance of each bin («{ ,f ). The 
X 2 value was computed as 


X 


2 


m 


2 


nf s 




(4) 


Table 3: Number of stars between 1.4 and 2.0kpc, and the ex¬ 
pected value from the exponential when we do not take into ac¬ 
count the arm bins. Significance and amplitude of the arm are 
also indicated for each of the three MB samples. 



°bs 

11 arm 

fit 

'‘■arm 

Signif. 

A 

MB-Slin 

303 

228 

4.3 

0.14 

MB-S2in 

302 

227 

4.3 

0.14 

CS-MBin 

182 

142 

3.0 

0.12 


where was computed from the surface density obtained from 
the fitted expression iff and transformed to the number of stars 

per bin using: nj ; " = iff -Sk- (Fyj). The factor {F zi ) is the aver¬ 
age of Fyj for all the stars in the bin previously used to transform 
from observed n k to observed surface density Y(r k ) (see Eq.QJ. 
Then, considering the degrees of freedom for each histogram 
(i.e., number of bins minus two, since Iir and^- 2 have been esti¬ 
mated), the test statistics were obtained. The hypothesis that the 
data come from a pure exponential distribution can be rejected at 
a 4-5cr confidence level for SI, S2, and CS samples for the inner 
sky area, so they clearly do not fit with an exponential. Looking 
at the plots, the overdensity at 1.6 kpc is the main reason for 
the deviation. The location of the peak of the distribution varies 
slightly for different samples, being at 1.5 kpc for MB-Slin, and 
1.7 kpc for MB-S2in and CS-MB in. An error of 0.2 kpc (the bin 
width) has been adopted for this maximum overdensity. 

On the other hand, when we rejected the points between 1.4 
and 2.0 kpc (where the overdensity is detected) and repeat the 
fit, we obtained p-values from the^- 2 test of 0.44, 0.50, and 0.10 
for the SI, S2, and CS samples (those in Fig. [3- In other words, 
when we did not take into account the bins located close to the 
arm the distributions were fully compatible with an exponential 
fit. 

We estimated the significance of the peak as {n ob r s m - 
nf rm )j where n ob r s m is the number of stars in the bins close 

to the arm (i.e., between 1.4 and 2.0 kpc) and n{ l ’ m is the ex¬ 
pected number of stars from the exponential fit (where the arm 
bins were not used for the fit). The overdensity obtained for the 
MB-S1 and MB-S2 samples had a significance of 4.3<x. For CS- 
MB we obtained 3.0cr, a lower value due to the smaller number 
of stars included in the sample (see Tabled. 

The amplitude of the arm A was estimated from the ex¬ 
pression A = (n* - riar m ) j{n'f r s m + nl‘ r ' m ). The obtained values 
are presented in Table [3] For all samples we derive values in 
the range A=[0.12,0.14], well in agreement with recent deter¬ 
minations (see references of the den sity c ontrast of the Milky 
Way spiral arms discussed in lAntoia et al.ll2Qllh . However, one 
shoul d take into account the diffe rent methods used by these au¬ 
thors: |DrimmeT&]S2ergeU20W^_obtained ,4=0.14 from K band 
surface brightness, while Beniamin et al.l (120051) derived A =0.13 
from GLIMPSE data looking at the tangential points_of the in¬ 
ner spiral arm. In external galaxies. iRix & Za ritskvl (1 1995 ) gave 
values in the range 0.15< A <0.6, while iGrosbol et aTH20041) 
provided amplitudes for several galaxies, also from surface pho¬ 
tometry in the K band, reaching values up to A =0.5, although 
they claimed that the amplit ude varies as a function of radius. 
However, the same authors (iGrosbol & Patsisll2014 .1 suggested 
that amplitudes measured in the K band may be overestimated 
by a factor of 2 which bring the Perseus amplitude in agreement 
with measurements of external galaxies. 
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Dist(pc) 



Fig. 4: Radial variation of the stellar surface density for the MB (top) and EC (bottom) samples in blue. Vertical red lines show 
the 1.2 and 3kpc completeness limits. The exponential fit is plotted in magenta, with the Hr and L 0 parameters expressed in pc 
and ★/pc 2 , respectively. Blue dots joined with solid lines are the ones used for the fit, since it is the region where completeness is 
ensured. 




Dist(pc) 


Fig. 5: Radial variation of the stellar surface density for the MB samples in blue. Vertical red lines show the 1.2 and 3 kpc complete¬ 
ness limits. Green dots show the points used for the exponential fit, avoiding those around the overdensity location. The exponential 
fit is plotted in magenta, with the Iir and S Q parameters expressed in pc and ★/pc 2 , respectively. 


6. The Perseus arm dust layer 

Since both individual photometric distances and visual interstel¬ 
lar absorptions (Ay) are available for a large sample of stars, a 
detailed 3D extinction map was created to discuss and quantify 
the possible existence of a dust layer related to the Perseus arm. 
In Sect 16.11 we present the maps and discuss some of the ob¬ 
served features, whereas in Sect. l6.2l the existence of a change in 
the differential absorption ( dAy/dr ) is discussed in order to de¬ 


termine if there is a high density interstellar region linked to the 
Perseus arm and, more importantly, to understand if this region 
is in front of or behind the stellar component of this arm. 
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distance (pc) 


Fig. 6: Number of stars n * observed (in blue) and number of stars 
expected from the exponential fit n J k " developed in Fig.0(in ma¬ 
genta) as a function of distance for the MBSlin sample. 


6.1. Three-dimensional extinction map in the anticenter 
direction 


The MB-S1 sample was used to create the 3D extinction map 
presented in Fig. [7] The number of stars in this sample allowed 
us to reach distances up to 2.5-3 kpc. The grid inside the 3D cone 
was constructed taking 20pc steps in distance and 2.5arcmin 
steps in both galactic longitude and latitude. At each point of the 
grid, the absorption Ay was computed as 


2£i ex P (iS) 

2Zi ex p(i^) 


(5) 


where N is the total number of stars in the sample and Ar, is 
the distance between the ith star and the point of the grid. The 
cr value is the radius used for the Gaussian weight computed 
as cr - r ■ atan(a) with a = 0.25°. This method allowed us 
to obtain a grid with higher spatial resolution in nearby regions 
where we had more information, and less at further distances 
where information was poorer. We note that this method would 
introduce some bias at the edges of the grid. The lack of data 
outside the surveyed area can result in mean values of the Ay lrlb) 
that are slightly biased towards the values present in the inner 
parts. We estimate that this bias should be present only in areas 
at angular distances less than about 0.5° (2<x) from the survey’s 
edges. 

In Fig. [7] we plotted the Ay values of the grid points placed 
at four distances and computed following Eq.0 A first compar¬ 
ison can be done with the integrated extinction maps obtained 
bv lFroebrich et al.l (120071) . using 2MASS data (Fig.0 top), and 
those recently derived by the Planck collaboration (2014; Fig.[8j 
bottom). As expected the same general trends are observed in 
these maps clearly identifying the same very low and very high 
extinction areas. Whereas both maps in Fig. [8] indicate the to¬ 
tal dust integrated intensity towards a given line of sight, our 3D 
map is able to distinguish features at specific distances. Although 
our data avoids the presence of some of the biases thoroughly 
discussed in lFroebrich et ah! d2007t) . we note that it is not the goal 
of this paper to investigate the individual distribution of clouds 
or their sizes and mass distribution, which are discussed in that 


paper. In the Froebrich maps, for example, regions close to rich 
star clusters would be dominated by the colours of the cluster 
members; instead, the data presented here do not suffer from 
this bias. As an example of the future scientific exploitation of 
this 3D map we here comment on the distribution of interstellar 
medium around the supernova remnant (SNR) Simieis 147. This 
SNR has its geometrical center at (l,b)=(18 0.17, b=-1.81) with 
a diameter of about 180arcmin (jGreenll2009l) . Its sugge_sted dis¬ 
tance ranges from 0.6 to 1.9 kpc (see Table 1 f rom Di ncel et al.l 
120151) . For its related pulsar PSRJ0538 +2817 iNg et alJ d2007l) 
estimated a distance of 1.47^2™kpc and lFesen et al.l ( 19851) an 


absorption of Ay = 0.76 + 0.2 magnitudes. At the edge of our 
survey our 3D extinction map indicates a slightly higher absorp¬ 
tion at larger distances (at about 2 kpc) and a more prominent 
absorption at th e area surrounding the geometrical center pro¬ 
posed bv [Green! (2009). As a second example we discuss the fea¬ 
tures visible arou nd (£, b ) ~ (1 82.0,0.0) in the 1.5 kpc and 2 kpc 
maps. iRodrfguez et al.l d2006l) found a molecular cloud in the 
anticenter direction associated with the two IRAS sources IRAS 


05431+2629 at (t,b) = (182.1,-1.1) and IRAS 05490+2658 at 
(C, b ) = (1 82.4,+0.3). This last HII region was also found by 
iBlitz et al.l ([1982 1 at ({, b) = (182.36, +0.19) through CO obser¬ 
vations, and they give a distance of 2.1+0.7 kpc. We observed a 
weak but well-defined high mean density region at this distance. 
These are only examples of the potential interest of these maps, 
and it will be through future analysis that multivariate clustering 
techniques will be applied to our grid points to confirm and char¬ 
acterize the clumpy structures present in the cone defined by our 
anticenter surve y. 

Recently, [Chen et al . (12014!) obtained a 3D extinction map 
with a spatial resolution of 3-9 arcmin by using multiband pho¬ 
tometry. In this study individual distances are not available, so 
to break the degeneracy between the intrinsic stellar colours and 
the amounts of extinction the authors relied on the combina¬ 
tion of optical and near-infrared photometric data. Our maps 
presented in Fig.[7J although with less spatial resolution, are 
based on individual distances with accuracies at about 10-20% 
jMonguio et al.l[2014 ). We understand that both studies provide 
complementary information for future analysis. 

As a final example, in Fig.[9]we show the 2D extinction map 
projected on the cartesian coordinates (x, y) and (x,z) planes 
where the absorptions at different z and y were averaged. The 
;c-axis is positive toward the anticenter, y in the direction of 
galactic rotation, while z toward the north galactic pole. These 
maps allowed us to analyze the change of extinction with dis¬ 
tance. We observed that larger extinction regions were placed 
at lower galactic longitudes and below the Galactic plane. This 
large-scale dust distribution is discussed in the next section in 
the context of the presence of the Perseus spiral arm. 


6.2. Absorption distribution as a function of distance 

Here, to reinforce the detection of the stellar overdensity of the 
Perseus arm presented in Sect. [3 we present in Fig.[T0](top) the 
median distribution of the visual absorption as a function of dis¬ 
tance. We can clearly appreciate a change in the slope at the same 
position where we detected the stellar overdensity of Perseus, 
which is at about 1.6 kpc. This change is better identified by plot¬ 
ting its derivative as a function of distance ( dAy/dr ) (see Fig. m 
bottom). We observed that this derivative presents two different 
well-defined values, a high value of about 0.8-1.0mag/kpc in 
front of the stellar arm, and a significantly lower value of about 
0.2-0.4mag/kpc behind the arm. We would expect, for constant 
absorption, a flat distribution of the dAy/dr , while the presence 


8 



















































M. Monguio et al.: First detection of the field star overdensity in the Perseus arm 
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Fig. 7: 2D extinction maps in ((, b) at four different distances from the Sun, 1000, 1500, 2000, and 2500. Ay is color coded, with a 
shift of 0.2 m in the color scale between consecutive plots. 


of a cloud or dust layer would be translated into a bump in the 
dAy/dr distribution. We can see this change in slope in all the 
directions, suggesting the presence of a dust layer just before the 
position of the density peak associated with the Perseus arm. 


As described by [Roberts .: (119691) . the star formation induced 
by the shock in a stellar density wave scenario would produce 
an azimuthal age sequence across the spiral arms. However, as a 
first attempt, we analyze the dust distribution (i.e., interstellar ex¬ 
tinction) in this context. The presence of the dust lane at the inner 
side of the stellar overdensity indicates a compression or shock 
in the gas associated with the spiral arm (although the strict co in¬ 
ciden ce of a shock and a lane is not obligatory uGittins & Clark j 
120041) . From this we can use our first analysis of the 3D extinc¬ 
tion map to establish a limit on the co-rotation (CR) radius in the 
Milky Way Galactic disk (' Roberts! 119721) . A discussion on the 
optical trace rs of spiral wave resonan c es in external galaxies can 
be found in lElmegreen et al.l (1 19921) . iPuerari & Dottoril (1 19971) 
provided a good sketch of the Z-trailing spatial arms configu¬ 
ration expected for the Milky Way (see their Fig. 1). Figure ITOl 
(right) suggests that the dust layer is just in front of the Perseus 
arm so the CR radius of the spiral pattern (if a density wave 
scenario is assumed) is outside the location of Perseus in the 
anticenter direction: Rcr > 10.2 kpc (assuming R a = 8.5 kpc). 
This result is closelyjn agreement with the position of the CR 
radius suggested bv lAntoia et al.l (1201 ll) . who, from a kinematic 
analysis of the moving groups in the solar neighborhood, pro¬ 
posed an angular velocity of the perturbation (T2 ;) ) of about 16- 
20km-s~ 1 kpW 1 , so with a galactocentric co-rotation radius at 
about Rcr-1 l-14kpc. 


7. Summary and conclusions 

In IMonguio et al.l d2013l) and IMonguio et al.l (12014l) we pub¬ 
lished a deep Stromgren photometric survey and a new strat¬ 
egy for the derivation of stellar physical parameters for the 
young stellar population in the anticenter direction. Here we 
used these data to derive samples of intermediate young stars, 
complete within the distance range between 1.2 and 3.0 kpc 
from the Sun. We computed the surface density distribution 
in the anticenter direction by assuming a vertical density dis¬ 
tribution of the disk, taking into account the surveyed area 
and the warp, so carefully defining correction functions for 
all the observational and physical effects. The same strategy 
was use d to compute the surf a ce den sity at the Sun’s position 
through lHauck & Mermilliodl ( 1998) Stromgren data, obtain¬ 
ing 0.032+0.001-fr/pc 2 and 0.017+0.001 ★/pc 2 for stars in the 
M v =[-0.9,1.2] and M v =[-0.9,0.8] ranges, respectively. The ex¬ 
ponential functions fitted to the computed stellar surface density 
revealed a clear star overdensity at around 1.6+0.2 kpc that we 
associate with the presence of the Perseus spiral arm. The ob¬ 
tained distance is sli ghtly closer than the 2.0 kpc recently ob¬ 
tained bv lReid et al.l ( 20141) who traced the arm using star form¬ 
ing regions. The fits we performed avoiding the bins close to 
the arm location, allowed us to derive the radial scale length of 
the Galactic disk young population, obtaining h K =2.9+0.1 kpc 
for B4-A1 stars, and h R =3.5+0.5 kpc for B4-A0 stars. The over¬ 
density associated with the Perseus arm has been detected with a 
significance above 3<x in all cases, and above 4cr for the samples 
with larger statistics. Our results indicate that the star density 
contrast of the young population in Perseus is on the order of 
,4=0.12-0.14. This quantitative estimation of the spiral arm am¬ 
plitude in the anticenter -derived here for the first time- should 
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Fig. 8: Top: Extinction map from iFroebrich et al.l d2007l) 
based on 2MASS data. Bottom: thermal the dust radi¬ 
ance R map (or dust integrated intensity) from Planck data 
(iPlanck Collaboration et al.]l2014l) . 
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Fig. 9: Top: 2D extinction map in (x, y), where stars for all z were 
used for the average. Bottom: 2D extinction map in ( x,z ), where 
absorption for the stars at all y were averaged. A v is color coded. 


be understood as a contribution to the future derivation of the 
amplitude changes of the Milky Way spiral pattern as a function 
of galactocentric radius. So far, this radial change has only been 
clearly detected in external galaxies fe.c.. Grosbol et all2004li . 

Also very important, the distribution of visual interstellar ab¬ 
sorption as a function of distance reveals the presence of a dust 
layer in front of the Perseus arm, which suggests that it is placed 
inside the co-rotation radius of the Milky Way spiral pattern. Our 
results locate the co-rotation radius at least two kiloparsecs out¬ 
side the Sun’s position. This is the first time that the presence of 
the Perseus arm has been detected through individual star counts, 
and its presence is being supported by the dust distribution. 
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Fig. 10: Top: Running median of the visual absorption vs. heliocentric distance distribution obtained from the three different samples 
MB-S1 (left), MB-S2 (middle), and CS (right). The gray area indicates the error of the median computed as 1.253 cr/ VfV. Bottom: 
differential visual absorption showing the change in slope of the upper plots. Vertical dashed lines show the location of the Perseus 
arm at 1.6 kpc derived in Sect. [4] 
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